Pressure effects on neutral and charged excitons in self-assembled InGaAs/GaAs 

quantum dots 
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By combining an atomistic pseudopotential method with the configuration interaction approach, 
we predict the pressure dependence of the binding energies of neutral and charged excitons: X° 
(neutral monoexciton) , X~ and X + (charged trions), and XX (biexciton) in lens-shaped, self- 
assembled Ino.6Gao.4As/GaAs quantum dots. We predict that (i) with applied pressure the binding 
energy of X° and X+ increases and that of X decreases, whereas the binding energy of XX 
is nearly pressure independent, (ii) Correlations have a small effect in the binding energy of X°, 
whereas they largely determine the binding energy of X~ , X + and XX . (iii) Correlations de- 
pend weakly on pressure; thus, the pressure dependence of the binding energies can be understood 
within the Hartree-Fock approximation and it is controlled by the pressure dependence of the direct 
Coulomb integrals J. Our results in (i) can thus be explained by noting that holes are more localized 
than electrons, so the Coulomb energies obey J (hh) > J (eh) > J (ee) . 
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The energetics of excitons reflects a balance between 
single-particle energy levels 8^ and of electrons 
(e) and holes (h) in the system, and the many-particle 
carrier-carrier interactions, resulting from electron-hole 
Coulomb and exchange interactionspi*^ The variation 
of excitonic energies under pressure naturally reflects 
the corresponding variations in single- vs many-particle 
energies. Of particular interest are the pressure varia- 
tions of excitons confined to nanosize dimensions such 
as in quantum dotsj 4 i 5 i 6 i 7 i 8 T 9 T 10 i 11 T 12 i 13 i 14 Unlike the case 
of excitons in higher-dimensional systems, where bind- 
ing and its pressure dependence reflects mostly many- 
particle (correlation) effects, in zero-dimensional (0D) 
systems where the geometric dimensions are smaller than 
the excitonic radius, binding of neutral and charged exci- 
tons results from an interesting interplay between single- 
particle and many-particle effects. Here, we use a realis- 
tic description of both single-particle and many-body ef- 
fects in self-assembled Ino.6Gao.4As/GaAs quantum dots, 
showing how pressure affects the different components of 
exciton binding. We distinguish the neutral monoexci- 
ton X° (one e, one h), from the neutral biexciton XX° 
(two e, two h), positive trion X + (one e, two h) and 
negative trion X~ (two e, one h). While the effect of 
pressure on X° has been measuredj&2ii2iiiii2iiii to the 
best of our knowledge, the optical spectroscopy oi X~, 
X + and XX under pressure has not yet been reported. 
For these reason, we provide definite predictions of the 
pressure effects. Each of the g-charged excitons has a 
spectrum of levels {v}, of which the lowest is termed 
the "ground state of x q " (x — X, XX). This spectrum 



is usually expressed by expanding the many-body exci- 
tonic states Y$^ v \x q )) y i a a se t °f Slater determinants 
|$(x 9 ))- The latter are constructed from single-particle 
electron and hole states and accommodate as many car- 
riers as are present in x q ■ The single-particle states are 
solutions to the effective Schrodinger equation 



1 , 
-2 V 



+ V ext (R) + V scr {R) \i>i = £ii>i, (f ) 



where V ext (R.) is the external (pseudo) potential (due to 
the ion- ion or ion-electron interaction) and V scr (R) is the 
screening response to such external potentials. The effect 
of pressure or strain is encoded in the ion-ion geometry 
underlying V ex t(R). The many-particle Hamiltonian is 
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where c\ (c,) and hj (hj) create (destroy) an electron 

in the single-particle state ip^ and a hole in ipj h \ re- 
spectively. In Eq. the Coulomb and electron-hole 
exchange matrix elements are given, respectively, by 
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Here, e(R, R') is a phenomenological, microscopic dielec- 
tric function that screens the Coulomb and exchange in- 
teractions, and in this work we have adopted the e(R, R') 

proposed by Resta^ 5 - 
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The diagonal elements ' 



; . , ( of Eq. J2J are the familiar electron-electron (/i/i' — 
ee), hole- hole (/i/i' = hh) and electron- hole (/i/i' = eft,) 
direct Coulomb integrals. The electron-electron and hole- 
hole exchange integrals are given by Jy.y and Jyly , re- 
spectively. Solving the single-particle Eq. Q for a given 
dot yields the wavefunctions ?/>i, which are used to con- 
struct the Slater determinants ^(x 9 )) f° r X q an d solve 
the many-particle, configuration interaction— (CI) prob- 
lem [Eq. (0)]. This gives the total (ground-state) energy 
Eci(x q ) °f exciton x q , as well as excitonic excited states. 

The binding energy of the excitonic complexes are de- 
fined as 
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In a simplified Hartree-Fock approximation and neglect- 
ing the electron-hole exchange ^oo'oo (whose magnitude 
is of the order of a few to hundreds of /j,eV)^ we have 
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The latter relation establishes a "sum rule" for the bind- 
ing energy of the biexciton at the Hartree-Fock level. The 
many-body, correlation effects 5(x q ) in the binding en- 
ergy can be quantified by comparing the full solutions in 
Eq. © to the HF ones in Eq. ©, 
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In this work, we consider a lens-shaped (base diame- 
ter b = 252 A and height h = 35 A) In . 6 Ga .4As/GaAs 
quantum dot and study how the excitonic binding ener- 
gies Aci(x q ) depend on pressure (well below the Ti c — 
Xq c crossover We then analyze this dependence in 
terms of the pressure dependence of (i) Coulomb inte- 
grals Jqq M ' and (ii) correlation energies 8{x q )- In Eq. 



we use a screened pseudopotential expressed as a 
superposition of screened atomic pseudopotentials 



^(R) + v scr (R) = v so + J2J2 < R - R i°"'> Tr ©]> 
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where Vso is a non-local spin-orbit pseudopotential; 19 
v a is a screened pseudopotential for atom of type a that 
depends on strain; and R; ' is the vector position of 
atom n of type a after the atomic positions within the 
simulation supercell (quantum dot+GaAs-matrix) have 
been relaxed, using a valence force field^S in order to 
minimize the elastic energy of the nanostructure. The 
explicit dependence of v a on strain transfers to the elec- 
tronic Hamiltonian the information on atomic displace- 
ments. v a has been fitted to bulk properties of GaAs 
and InAs, including bulk band structures, experimental 
deformation potentials and effective masses, as well as 
LDA-determined band offsets^ Equation is solved in 

a basis of linear combination of Bloch bands {w^jf (R> ^)} 
with band index n and wave vector k of material M 
(=GaAs, InAs) strained^ , to e. Thus, 
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where N is the number of primary cells in the simulation 
supercell that contains the quantum dot and GaAs ma- 
trix. The many-body configuration-interaction (CI) ex- 
pansion is taken over all the Slater determinants \&(x q )) 
generated within a set of 12 electron and 20 hole single- 
particle, confined states. 

Direct- Coulomb vs correlation contributions to 
binding. — Table [I] shows the Cl-calculated binding 
energies Aci{x q ) as well as its decomposition [Eq. 
0] into Hartree-Fock AnF{x q ) and correlation 5(x q ) 
contributions. We see that (i) the binding energy of the 
neutral monoexciton X° is constituted primarily by HF 
energy with only 6 % being due to correlation. This is 
in contrast with X a in bulk semiconductors where S(X°) 
dominates over A HF (X°). (ii) For X~, X+ and XX 
the HF and correlation contributions to binding are 
comparable. Specifically, while X~ is bound (positive 
A) already in HF, here X + and XX are unbound in 
HF, but become bound by correlation, (iii) For each 
excitonic complex, the magnitude of the correlations 
depends weakly on pressure. 
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TABLE I: Comparison of Hartree-Fock (HF) and many-body configuration-interaction (CI) binding energies (in meV), and 
verification of the "sum rule" [Ahf(X~) + Ahf(X + )] for different pressures. For each excitonic complex x Q > we present 
the CI binding energy Aci{x 9 ) as a sum of the Hartree-Fock binding energy AHF(x q ) an d correlation-energy component 
<5(x 9 )=Ac/(x , )-Ahf(x 9 )- 





Quantity 


0.2 GPa 


0.8 GPa 


1.3 GPa 


1.8 GPa 


2.4 GPa 


A C /(X U ) 


= A HF (X»)+S(X U ) 


20.8 + 1.3 


21.1 + 1.4 


21.3 - 


hl.4 


21.6 + 1.5 


21.9 + 1.5 


Aci(X-) 


= Ahf(X-)+S(X-) 


1.2 + 1.3 


1.0 + 1.4 


0.8- 


-1.5 


0.6 + 1.6 


0.5 + 1.6 


A C i(X+) 


= A hf (X+) + S(X + ) 


-1.8 + 2.4 


-1.5 + 2.3 


-1.2 - 


-2.3 


-0.9 + 2.2 


-0.7 + 2.2 


A c i (XX 


) ) = A H f(XX°) + 5(XX°) 


-0.6 + 2.0 


-0.4 + 2.0 


-0.3- 


-2.0 


-0.3 + 2.0 


-0.2 + 2.0 


Aci(X-) 


+ A C i(X + ) 


3.1 


3.2 




3.4 


3.5 


3.6 


Sum rule 




-0.6 


-0.5 




-0.4 


-0.3 


-0.2 



Pressure dependence. — Figure^shows the dependence 
on pressure of (a) the binding energies Ac/(x 9 ), (b) 
Coulomb energies Jqq 6 ' , Jqq and Jqq , and (c) correla- 
tion energies S(x q )- Pressure is represented by Aa/ao = 
(a — arj)/an, where a and an are the distorted and 
equilibrium lattice parameter of the GaAs matrix, re- 
spectively. The pressure values showed in the upper 
axis on Fig. Qfa) are calculated by using the equa- 
tion of state 21 P = (Bo/B' )[(V /V) B '° - 1], where we 
take V /V = [1 + Tr(e)]- 1 and calculate Tr(e) in the 
GaAs matrix away from the dot. The calculation of 
e = e(R) is performed using atomistic elasticity^ We 
take B = 74.7 GPa and B' Q = 4.67 as the GaAs bulk 
modulus and its derivative with respect to pressure, 
respectively^ We see from Fig. Q] that the pressure de- 
pendence of the binding energy of the various excitons 
is different: (i) Aci(X°) shows a small, nearly linear 
increase with pressure; changing about 7 % in the stud- 
ied pressure range. Aa(X~) decreases slightly with in- 
creasing pressure, while Aqi(X + ) increases significantly; 
at Aa/a = -0.0087 (P = 2.4 GPa) it has increased by 
160% compared to its value at Aa/ao = 0. Similar to 
the monoexciton case, the binding energy of the biexci- 
ton depends only weakly on pressure, showing a small 
relative change as pressure reaches 2.4 GPa. (ii) Equa- 
tion J7J) shows that the binding has a HF part and a 
correlation part. Table showed that the magnitude of 
the binding is decided by the HF part for X° and by 
both HF and correlation for X~ , X + and XX . How- 
ever, Fig. ^shows that the pressure dependence is always 
decided by the HF contribution, (iii) The "sum rule" of 
Eq. |[BJ valid within HF is not valid at the CI level, 
quantitatively failing to predict the correct values of the 
biexciton binding energies, (iv) From (ii), we see that 
the trends of the binding energies with pressure are de- 
termined by Jqq , JqI^ and Jqq ^ . By calculating these 



integrals, we find that Jqq^ 



> J, 



(eh) 



> Jqq^ and the mag- 



nitude of Jqq and J^f 1 increase with a similar slope as 
pressure increases while Jqq remains nearly constant; see 
Fig. This explains the decrease of Aci(X~) and the 
increase of Aci(X + ) with applied pressure. Further, it 
also becomes clear why the binding energy of the biex- 
citon remains nearly unchanged with changing pressure: 
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FIG. 1: Pressure dependence of (a) binding energies Aci(x 9 ) 
as obtained from the many-particle, configuration-interaction 
method, (b) Coulomb-energy component Jqq'' ' and (c) 
correlation-energy components S(x q ). (See text for defini- 
tions.) 



The similar rate of increase of Jqq 6 '' and Jqq 1 ' with pres- 
sure combined with the magnitude of Jqq leads to a 
weakly pressure-dependent binding energy for XX . 

Wavefunction localization with pressure. — To under- 



T (eh) 



stand the trend Jqq ^ > J™ > J 



< - ee - ) and also that Jqq ^ 



oo 



has the weakest pressure dependence, while Jqo^ has the 
strongest, Figure [21 shows the calculated wavefunctions 
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LUMO and HOMO wavefunctions vs pressure 
Lens-shaped lr\, 6 Ga 04 As/GaAs; b=252A, h=35A 




FIG. 2: LUMO [V>o e) ] and HOMO [ip ( h) ] wavefunctions at 
0.2, 1.3 and 2.4 GPa. The outline of the dot is present as a 
light shadow. The (red) isosurface encloses 75% of the charge, 
the in-plane (bottom) contour plot is taken at 1 nm above the 
dot's base and the out-of-plane contour plot bisects the dot. 
The LUMO wavefunction becomes more confined as pressure 
increases — see the how the wavefunction penetrates less into 
the barrier above of the dot as pressure increases; whereas the 
HOMO wavefunctions is almost independent of pressure. 



for the electron ground state [ipQ^] and hole ground state 

[tpQ 1 ^] as a function of pressure. The isosurfaces enclose 
75% of the charge density, the in-plane contour plot is 
taken at 1 nm above the base of the dot and the out- 
of-planc contour plot bisects the dot. We see that the 
electron is always more localized than the hole. In addi- 
tion, the electron gets more localized as pressure is ap- 
plied, while the localization of the hole remains nearly 
unchanged. The in-plane (parallel to the base) spatial 
cxtcntion of the electron does not change as much as the 
out-of-plane. In particular, Fig. [21 clearly shows that 
the spatial penetration of the electron wavefunction into 
the GaAs matrix decreases with applied pressure. The 
increased localization of the electron with pressure can 
be explained by the larger magnitude of the conduction- 
band edge (CBM) deformation potential of bulk GaAs 
with respect to that of bulk InAs, 24 which results in an 
increased electron confinement in the dot with pressure. 
In contrast, the similar magnitude of the valence-band 
edge (VBM) deformation potential of both bulk GaAs 



and InAs leads to small changes in hole confinement with 
pressure and, therefore, to small changes in localization. 

Single-particle and excitonic pressure coefficients. — 
We calculate the linear pressure coefficient a by fitting 
the pressure dependence of the band gap to E g (P) = 
E g (0)+aP+bP 2 . For the dot, at the single-particle (SP) 



level, we obtain 
excitonic value is 



(dot) _ 



SP 



(dot), 
l CI 



86.47 (meV/GPa), whereas the 



(A -0 ) = 85.79 (meV/GPa). The 
latter compares well with the values of 85 (meV/GPa) 8 , 
80(meV/GPa)ii and 82 (meV/GPa) 11 observed in 
InAs/GaAs dots for the emission lines at 1.28eV, 1.26eV 
and 1.30 eV, respectively. For bulk GaAs, we obtain 
a (buik) _ 205.86 (meV/GPa), which is within the range of 
observed values: 94-120 (meV/GPa)i2i2i Thus, the dot 
has a smaller linear pressure coefficient than bulk GaAs. 
In addition, it is interesting to inspect how the VBM 
and CBM contribute to the linear pressure coefficient of 
the band gap. By fitting the lattice-deformation (pres- 
sure) dependence of single-particle eigenvalues £j(a) = 



(dot) 
VBM 



£j(ao) + Ai(Aa/a ) + B^Aa/ao) 2 , we find A 
-1.76 eV and = -22.49 eV. For, bulk GaAs we 

find A vbm = -3.74 eV and ^ CBM = -27.23 eV. We 
see that the band-gap response to the lattice distortion 
(pressure) is largely dominated by the changes in CBM. 
To reproduce this LDA-predicted 2 ^ behavior, it is neces- 
sary to have a pseudopotential that explicitly depends 
on strain, otherwise one gets ^4cbm = — 0.46eV and 
^4vbm = — 11.72 eV and, consequently, VBM dominates 
the gap changes. 

In summary, we have studied the effects of pressure 
on the binding energies of X°, X~ , X + and XX . Our 
main findings are the following, (i) With applied pres- 
sure, the binding energy of X° and X + increases and that 
of X~ decreases, whereas the binding energy of XX is 
nearly pressure independent, (ii) The correlation-energy 
component in the binding energy of X° is small, whereas 
it is large in X~, X + and XX ; indeed, correlation is 
fully responsible for binding the latter complexes, (iii) 
Correlations depend weakly on pressure, (iv) The pres- 
sure dependence of the binding energies is controlled by 
the pressure dependence of the direct Coulomb integrals. 
Further, the relative magnitude (order) of these direct 
integrals explains the relative magnitude (order) of the 
binding energies, (v) Pressure dependence of J^ h \ Jqq^ 

and Jot) 6 '* is explained by the changes of the LUMO and 
HOMO wavefunctions with pressure. 
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